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Abstract 

We derive a simple formula for the stopping distance for a high-energy quark traveling through a 
weakly-coupled quark gluon plasma. The result is given to next-to-leading-order in an expansion in 
inverse logarithms ln(£'/T), where T is the temperature of the plasma. We also define a stopping 
distance for gluons and give a leading- log result. Discussion of stopping distance has a theoretical 
advantage over discussion of energy loss rates in that stopping distances can be generalized to the 
case of strong coupling, where one may not speak of individual partons. 
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I. INTRODUCTION AND RESULTS 



There has long been interest in calculating how a fast moving quark or gluon looses 
energy when traveling through a quark-gluon plasma. One of the simpler versions of this 
problem is to calculate energy loss in the case of a spatially-infinite equilibrium plasma (i.e. a 
plasma whose size is large compared to the stopping distance). For weakly-coupled plasmas, 
energy loss is dominated by hard bremsstrahlung and pair production in the high-energy 
limit E ^ T. Results for the rates of these processes at leading order in have been known 
for some time, as we'll review in a moment. The rate is parametrically 0, 0, 0] 
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for hard bremsstrahlung and pair production at high energy (with some caveats on the 
range of applicability to be discussed later). The roughly E'^^"^ decrease of this rate with 
energy is a result of the Landau-Pomeranchuk-Migdal (LPM) effect [gI, fk^ , which arises from 
interference between bremsstrahlung (or pair production) associated with multiple collisions 
of the particle with the plasma. If there were no LPM effect, then the high-energy rate 
would have been approximately independent of E. 

From the rate of bremsstrahlung and pair production, one can figure out how far a high- 
energy particle travels before it "stops" in the plasma, by which we mean that the particle 
loses enough energy to reach ~ T and so equilibrate with the plasma. Consider the 
energy loss of a quark. Crudely, if the quark loses roughly half its energy in the first hard 
bremsstrahlung, and roughly half in the next, and so forth, then the roughly E"^^"^ behavior 
of the rate (11.11) suggests that the stopping distance is parametrically of order 
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In weak coupling, the E~^^'^ fall-off of the rate is due to the LPM effect. It's interesting to 
ascertain whether anything similar happens in plasmas that are not weakly coupled. Fortu- 
nately, there are certain strongly-coupled gauge theories where people have, using AdS/CFT 
duality, calculated results related to the slowing down of high-momentum particles. We can- 
not directly discuss bremsstrahlung or pair production rates at strong couphng, however, 
because bremsstrahlung and pair production implicitly refer to individual quanta and so are 
intrinsically perturbative concepts. But one can generalize the idea of stopping distances 
to non-perturbative situations, as nicely explained by Chesler, Jensen, Karch, and Yaffe jl] , 
as we shall review later. In the strongly-coupled gauge theories that have be en st udied, 
stopping distances at high energy behave like E^^^ ]1, M, 13] rather than E^^^/ Vln E. 



Because stopping distances, unlike bremsstrahlung rates, can be generalized to strongly- 
coupled situations, it's interesting to know how to calculate them in weakly-coupled situa- 
tions. In this paper, we derive simple results for various stopping distances in QCD in the 
high-energy limit ln(£^/T) ^ 1, working to leading order in coupling ctg. Though we assume 
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TABLE I: Coefficients a,b,c for the quark-number stopping distance p.3|) in the case Q± ^ T 
(parametricahy E ^ T^/q or equivalently E ^ T/a'^ln{a^^)). Also shown are the leading-log 
coefficients Og'^'' and Oq"'' for the gluon and quark energy stopping distances defined in Sec. IIIB[ 
Note that larger values of a correspond to shorter stopping distances in the high-energy limit. 
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TABLE IL Coefficients for the stopping distance in the case Q± <^ T (parametricahy T <^ E <^ 
T^/q^T/alHa-^)). 



\a{E/T) ^ 1, we will formally assume that as is small enough that Q;sln(£'/T) <^ 1 and so, 
for instance, we ignore the running of the coupling a^- (We will say a few words later about 
how to plausibly modify the answer to include running coupling effects.) 

Our most thorough calculation is of the stopping distance for a quark, which we compute 
to next-to-leading-logarithmic order (NLLO) in \\i{E/T). Our results can be summarized 
as 

1 

4top,q^ 2^V^' 

aai_i \ 1 L 

where the logarithm L is determined by self-consistent solution of the equation 

and a, 6, and c are numerical constants. These constants are given in Table [I] in the large 
E limit of E > T/al\Yi{a-^), or in Table [III in the case T <^ E <^T /a^ hii^a^ ^) of smaller 
(but still formally large) E. The latter case has c = in (]1.3bp and so corresponds to the 
simple parametric estimate fll.2p based on fll.ip . 

Later, we will also give a definition of a stopping distance for gluons that can be gener- 
alized to non-perturbative gauge theories, and we will give a result to leading-log order for 
that distance in the case of weak coupling. 
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At the moment, we do not have a precise weakly-coupled result for the types of super- 
symmetric gauge theories in which the strongly-coupled limit has been studied. The answer 
will be in the form (I1.3P for weak coupling, but we leave the calculation of the coefficients 
a, b, and c in this case for future work. 

Throughout this paper, we will assume that the baryon chemical potential of the quark- 
gluon plasma is negligible compared to its temperature. 

In section ini we review Chesler et al.'s definition of a stopping distance for quarks, which 
we call the quark number stopping distance, and then we discuss a natural generalization for 
discussing gluons, which introduces what we will call the quark and gluon energy stopping 
distances. Section [TTl] derives results for these various stopping distances to leading-log 
order, which is enough to determine the coefficients a and c in our resuh ([L3D. Section HV] 
then continues to next-to-leading log order for the quark number stopping distance, giving 
the remaining coefficient b. As a test of our NLLO result, we compare it in section |V] 
to numerical simulations of energy loss that make no approximation regarding the size of 
ln{E/T) but instead use full weak-coupling results for the bremsstrahlung rate. We verify 
that our NLLO result (11. 3p correctly reproduces the large E behavior. Finally, in section 
IVit we discuss how one might modify our result to account for running of Og. 

II. DEFINING STOPPING DISTANCES 
A. Quark Number Stopping 

To explain Chesler et al.'s definition of a stopping distance for quarks, we first focus on 
the case of weak coupling, where we can talk about individual particles. We gave a crude 
estimate of the stopping distance in flL2l) . but we should consider that an individual quark 
doesn't really stop in a weakly-coupled plasma. Once its energy drops to -E ~ T, then the 
quark is just as likely to gain energy from its interactions with the plasma as lose energy. 
A cartoon of an individual high-energy quark's trajectory is shown in Fig. [T^. At first, 
the quark moves in a nearly straight line as it loses energy. But once its energy falls to 
be ~ T, it then random walks through the plasma just like any other equilibrated quark. 
Suppose we repeat this thought experiment over and over again in order to determine the 
probability distribution of the quark's position as a function of time. A cartoon of the time 
evolution of this distribution is shown in Fig. [T)d. Once the quark equilibrates, the center of 
the distribution stops moving, and the distribution merely spreads in size due to diffusion. 
We can then define the stopping distance of a quark as the distance that the center of the 
probability distribution moves between t=0 and t=oo. 

So far, we have been speaking perturbatively about the position of the quark. The 
probability distribution of quark positions easily generalizes to the non-perturbative case: 
simply refer instead to the probability distribution for total quark number in the system, 
since the quark number operator is defined non-perturbatively. But so far we've specified our 
initial state perturbatively when we said it was a single high-energy quark moving through 
the quark-gluon plasma. Chesler et al. realized that this can also be generalized to the non- 
perturbative case. Consider all possible initial states that (i) differ from thermal equilibrium 
in a localized region of space, (ii) contain total energy E in that region in excess of the 
thermal background, and (iii) contain total quark number one in that region. Measure the 
stopping distance for that initial state. Then define the "quark number stopping distance" 
as the maximum stopping distance, taken over all such choices of initial state. 
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stopping distance 



FIG. 1: The stopping of a high-energy quark represented by (a) the path of an individual quark, 
and (b) the motion of the center of the spreading quark number probabihty distribution. In (b), the 
quark number probabihty distribution at each successive time is represented by a single contour. 
As time progresses, the center of the contour moves from left to right, slowing down with time, 
and the diameter of the contour increases. 

Let's understand how this definition would work out in the weakly-coupled case. In our 
earlier discussion, our initial state was a single quark carrying the total energy E. What if 
we had chosen a different initial state with the same energy and quark number? Consider, 
for instance, an initial state composed of a quark of energy E/2 plus a gluon of energy E/2. 
The quark number stopping distance doesn't care about what happens to the gluon and 
would just depend on how long it took the energy E/2 quark to stop. But according to 
f ll.2p . that is roughly l/\/2 shorter time than if the entire energy E had been carried by a 
single quark. So the initial state that gives the maximum stopping distance is indeed the 
one we want in the weakly-coupled case: all the energy carried by a single quark. 

The definition of stopping time just given implicitly assumed that both the position and 
energy of the initial state are well defined. Due to the uncertainty principle, however, the 
location of the initial position will naturally be ambiguous by an amount of order 1/E. So 
the definition of stopping distance implicitly assumes the limit that £stop ^ l/E. 

In the case of A/'=4 large- A^^c supersymmetric Yang Mills theory with a massless Af = 2 
fundamental-charge "quark" hypermultiplet, Chesler et al. find (with their non-perturbative 
definition of the stopping distance) 



B. Gluon Energy Stopping 

Even in weak coupling, there is some difficulty in defining exactly what one means by 
the stopping distance of a gluon. For quarks, it's always clear what the energy of the quark 




(2.1) 



in the A = Ncg^ 



oo limit, with 



C ^ 0.5. 



(2.2) 
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is after a. q ^ qg splitting process. For gluons, we can have g ^ gg or g qq, and we 
have to decide which particle to follow. One possibility would be to follow, at each splitting, 
the final particle that has the highest energy. This definition does not have any obvious 
generalization to strong coupling. But we can instead use a definition like the previous one 
for quark number stopping distance. For the initial state, we now require quark number 
zero. Then, for the stopping distance, we measure the distance traveled before the center of 
the probability distribution for energy in excess of equilibrium (rather than quark number) 
stops. We will call this the "gluon energy stopping distance." 

For a gluon in a weakly-coupled plasma, this definition can be translated as follows. 
Consider the splitting of an initial high-energy gluon, moving in the z direction, which 
cascades through splitting into particles, which are "stopped" at positions zi, Z2, zn 
relative to the initial position of the initial gluon. Here stopped means that the individual 
energies Ei are of order T. The gluon energy stopping distance is the energy-weighted 
average of these stopping positions: 



^(energy) ^ Yli ^i^i 
f-stop,g — 



(2.3) 



This simple translation is fuzzy when it comes to deciding exactly when a final parton has 
low enough energy to be considered "stopped." Is it £'j < T or i?j < 2T or ...? However, 
because the time scale (a^T)~^ for particles with energy of order T to significantly interact 
is so much faster than the time scale (II. 2p for the initial gluon to lose significant energy, 
this ambiguity only involves corrections to the calculation of the stopping distance that are 
suppressed by a factor of \JT/ E (up to logarithms). We can ignore these corrections if we 
are just interested m E ^ T . Note that, when E^ ~ T, we should also be thinking about 
collisional energy loss, which is then competitive with bremsstrahlung and pair production. 
But the effects of collisional loss on the stopping distance will similarly be suppressed by a 
factor of \JT/E. 

III. LEADING LOG 

A. Quark Number Stopping: General Analysis 

We'll begin with a leading-log calculation of the quark number stopping distance, which 
is fairly simple. To motivate the approximation, consider the crude, heuristic estimate made 
in (11. 2p . At leading-log order, we can ignore the difference between a hi{E/T) appearing 
in T{E), a \n{E/2T) appearing in r(ii^/2), and a hi^E/AT) appearing in r(i?/4). Since the 
series in (II. 2p converges sufficiently rapidly (i.e. geometrically), it is sufficient at leading-log 
order to replace the logarithm in the bremsstrahlung rate (II. ip by a constant: 



nE)^ciT^'l^ (3.1) 

where Eq is the initial energy of the quark. In this approximation, the E dependence of 
T{E) is simply E~^^'^, without any logarithm. That simplification will be the key to a simple 
solution for the stopping distance in what follows. 

In more detail, the bremsstrahlung rate is a function of the momentum fraction x = u/E 
of the emitted gluon of frequency u. (For hard bremsstrahlung in the large E limit, we can 
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ignore the thermal mass of the gluon.) In the leading- log approximation, we can write 

dx \ E J dx dx ^ 

where there is an implicit logarithmic dependence of dV/dx on Eq. We will review the 
specific leading-log formula for q ^ gq bremsstrahlung later. But for now we can proceed 
with the analysis in a general way, based simply on the scaling (13. 2p . 

Let iq{E) be the quark number stopping distance as a function of E. We can formally 
write the following self-consistent equation for i^. 



.(^) = T^-^+ fdx ^^7^^^;;|/^" f,((l-a:)E). 



(3.3) 



The first term on the right-hand side is the average distance before the first bremsstrahlung. 
That first bremsstrahlung splits the quark into a nearly-collinear gluon with energy xE 
and quark with energy (1 — x)E. The second term on the right-hand side of (13. 3p is the 
remainder of the stopping distance after the first bremsstrahlung. This is just ^q((l — x)E) 
weighted by the probability [dV /dx)/T that the first bremsstrahlung had gluon momentum 
fraction x. To find the stopping distance, we just need to solve this equation in leading-log 
approximation. 

Note that (13.31) ignores what happens to the radiated gluons. Suppose a bremsstrahlung 
gluon later splits into a qq pair with the new quark carrying more energy than the new 
anti-quark. The subsequent cascading of that pair will then affect the distribution of quark 
charge at the end of the cascade process, and so will affect where the center of the total quark 
number distribution stops. However, it's just as likely that the gluon would have instead 
split into a qq pair in which the new anti-quark carries more energy, and this will produce an 
opposite effect on the location of the center of the final quark number distribution. Because 
of charge symmetry, we can ignore what happens to the gluons when calculating the quark 
number stopping distance. We dwell on this point because the situation will be different 
when we later compute the gluon energy stopping distance. 

The reason we referred to (13. 3p as "formal" is that the total bremsstrahlung rate 

Wi?)=^'rfx^^^^^=^^ (3.4) 

is infrared (x— >-0) divergent in leading-log approximation. You can imagine temporarily 
imposing some sort of infrared regulator. But it is easy to manipulate (13. 3p to remove the 
need for such a regulator by multiplying both sides by Vq^gq{E) and rearranging terms to 
rewrite the equation as 

dTq^gq{E,x) ^^^^^^ _ ^^^^^ _ ^^^^^ ^ ^ ^^^^^ 

The fact that the factor in brackets vanishes as x — > turns out to be sufficient to make the 
x integration in (13.50 infrared convergent. 

Now make the leading-log approximation by substituting (13.21) into (13. 5p to get 

E''/' dx ^^^^^^^ [UE) - U{1 - x)E)] = 1. (3.6) 



7 



We can solve this equation by writing 

^1/2 



'c(^) = (3.7) 



where is a constant. Then the factors of E drop out of (13.61) . leaving 

= f dx [1 - (1 - xy/^] . (3.8) 

Jo dx 

It only remains to plug in the detailed leading-log formula for dTg^gg/dx in leading-log 
approximation. We'll leave that to Sec. IIII C[ 



B. Gluon Energy Stopping: General Analysis 



For the gluon energy stopping distance ig (E) given in weak coupling (and high energy) 
by (12. 3p . we can write down an equation similar to (13. 3p . One difference is that we now have 
to account for both of the high-energy particles present after the first splitting, because they 
both carry energy. To warm up, first consider pure Yang Mills gauge theory with no quarks. 
Then the only relevant splitting process would he g gg, and the equation for the gluon 
energy stopping distance would be 

^t\E) = -^ + l j\ /^9^9,^E x)/dx [^^(e)(^^)^(i_^)^(e)((i_^)^)] (3.9) 

The new feature to this equation is that the energy stopping distances ^g\xE) and ^g'^ ((1 — 
x)E^ of the two gluons after the first splitting are averaged, weighted by their energies, as 
required by (12.31) . (An energy- weighted average of the energy stopping distances of these 
two gluons is the same thing as an energy-weighted average of the positions of all their 
stopped descendants.) The overall factor of | in front of the integral in (13.91) is to avoid 
double counting states of the final two, identical gluons. 

Once quarks are introduced, there is the added complication that the energy carried by 
a gluon could be converted into energy carried by quarks via g qq. To find the gluon 
energy stopping distance i^g\E), we will therefore need to also introduce the quark energy 
stopping distance iq\E), which is given in weak coupling by ( 12. 3p but in the case where the 
initial particle is a quark (or anti-quark). The basic equation (13. 9p for i'^^ now becomes 

^i\E) = Y^^+iy^{l "^^^-f^^^^^/"^' [xi^ixE) + {l-x) ((1 - x)E)] 

+ ^^^-f^^'^''^'' \x^{xE) + (1 - X) £f) ((1 - x)E)]Y (3.10) 

where 

r, - r,^,, + r,^,, = \j\x ^ + dx ^ (3.11) 
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is the total rate for a gluon to split by either bremsstrahlung or pair production. We now 
need the corresponding equation for the quark energy stopping distance: 

To solve this coupled system of equations, we proceed as before. First eliminate the issue 
of infrared divergences by rewriting them as 

' 1 [£(e)(^) -xi^^\xE) - (1 - x) ((1 - x)E)] 

+ '^^'^'^^^'"'^ -xi^\xE) - (1 - x) £f) ((1 - x)E)] 1=1, (3.13) 



dx '^^ ''~'9<](^ ' ^ 



'^-[l^;i\E)-xl^^\xE)-{l-x)£^;i\{l-x)E)] =1. (3.14) 



dx 

Then we make the leading-log approximation fl3.2p and write 



4°H^) = ^' ^'^^E) = ^^ (3-15) 



q 

to obtain coupled algebraic equations for l/Ag'^ and 1/A^^\ The solution is 



where 



1 ^I^^mM [1 _ a;3/2 _ (1 _ ^)3/2] ^ ^fl^^^^l^ (3.16b) 

Mgq = /' dx '^l^d^ [_a;3/2 _ (1 _ a;)3/2] ^ (3_lg^) 







M,, = rfx ^5:^^^ [1 - (1 - :,f/^] . (3.16e) 



dx 

C. Leading Log Details 

1. Bremsstrahlung and pair production rates 

Following notation similar to Ref. the leading- log splitting rates are 



dVa^bc _ a iil^^^P^^^,{x) 



dx 4nV2x(l - x)E 



(3.17) 
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where, at leading-log order, 



|4x(l - x)E[{-Ca + a + C,) + {Ca -C, + C,)x^ 

+ (a + a - - x)2]t(gxo)}'^'. (s.is) 



Here ^(A) is proportional to the average squared transverse momentum Q\ that a high- 
energy particle picks up per unit length while traveling through the plasma. We'll discuss the 
formula shortly. In the present context, the definition of ^(A) is UV-regulated by imposing 
an ultra-violet cut-off A imposed on the momentum transfers from individual collisions. 
Q±o ~ {qEY^^ is any rough estimate of the total transverse momentum transfer during one 
bremsstrahlung formation time,^ and ambiguities in that guess of 0(1) factors only affect 
the answer beyond leading-log order. Cs is the quadratic Casimir associated with the color 
representation of particle s. For QCD, 

Cq = Cy = |, Cg = Ca = 3. (3.19) 

The Pa^b are the usual Dokshitzer-Gribov-Lipatov-Alterelli-Parisi (DGLAP) splitting func- 
tions 

Pq^g(x) = Cf ^^^^^7""^'^ , (3.20) 

Pg^q(a;) = mF[x^ + (1 - x)% (3.22) 

where iVf is the number of quark flavors and tg is the trace normalization of color generators 
T" deflned by tT{T^T^) = with 

= I , tA = Ca, (3.23) 
The UV-regulated ^(A) discussed above is given in weak coupling by 

g.(A) ^ C.t(A) = / d'qs. ^ ql, (3.24) 

Jq^<A d^(l± 

where Fei is the rate of elastic 2— >2 scattering of a high-energy particle of type s off of the 
plasma, and q± is the transverse momentum transfer in an individual such collision. The 
only dependence on the type of high-energy particle is a factor of Cg, which we have factored 
out of the deflnition of q. 



1 



Here is a lightning review. By the uncertainty principle, the formation time tf is of order l/AE, where AE 
is the amount by which energy would be violated if an on-shell high-energy particle split by bremsstrahlung 
(or pair production) in isolation. For hard bremsstrahlung, AE ^ Q\/ E in the high-energy limit, where 
Q±^ is the transverse momentum of the final particles, and so if ^ ElQ\. The must be supplied by 
collisions with the plasma. By definition of g, the amount of Q\_ picked up in one formation time is given 
by Q\ ^ <jtf- Putting the last two equations together gives t{ ~ {E/qY/'^ and ^ {qEY^^. (For a 



review, see, for example. Sec. 3.2 of Ref. [l3| together with Ref. [14|.) 
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The weak-coupling result for q{A) is slightly different depending on the precise range of 
A. At leading-log order in coupling,^ 



1(A) ~ aTml In (^^^ + 47raWln (^^^ 



if A > T, where the second term dominates in the limit of large A; and 

/ A^ 

q(A) ~ aTmi In 

if mj) ^ A < T. Here is the Debye mass, given by 

ml = i(tA + iVftp)/T2 = (1 + 1 N,)g'T\ 
and M is the plasma particle density weighted by group factors, given by 

AT = |(tA + I iVftp)|| = (1 + i 
where ({s) is the Riemann zeta function. 



(3.25) 



(3.26) 



(3.27) 



(3.28) 



2. Leading log stopping distances 



Substituting the leading-log bremsstrahlung and pair production rates into the general 
leading-log formulas (13. 8p and (13.161) gives us leading-log stopping distances. For the quark 
number stopping distance (13.81) . 



a 



Aq = — -=[|(g^o)]'/' / dxP^^,i 
Zn\/2 Jo 



X) 



Ca + (2Cf - Ca)x^ + Ca(1 - x) 
x{l — x) 



21 1/2 



Now take Qxo ~ (qE) 



1/4 



X [1 - (1 

■T^Ef'^. If Qxo < T, then fl3:26|) gives 



X 



q{Q 



aTm-Q In , „ 
mi 



~ ^aTml In 



at leading-log order. Note that the factors of coupling cancel out in the logarithm. If 
Q±o ^ T, and if we keep track of logarithms of coupling as well as logarithms of energy, 
then 



T 



], (3.29) 



(3.30) 



q{Q±o) - aTml In 



+ Area' AT In 
mn I \ 



In 



«-V|.„(il 

2 V"^ 



m 



D 



D / J 



ln(f)+cln(i 



27raWln 



E 



(3.31) 



^ For a discussion of the perturbative result for (j(A) using this notation, see Ref. 
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with 



c = 2 



Tml 



- 1 



In what foUows, it wiU be convenient to define the dimensionless numbers 



mo = 



aT2 



{tA + Nitp), 



We can now summarize the leading-log result as being given by f ll.Sp with 



(3.32) 

(3.33) 
(3.34) 



Air 



[x] 



Ca + {2Cf - Ca)x^ + Ca{1-x) 



Q = 



x(l — x) 

ml, < T, 



1/2 



and 



2(^-1), g±»T; 



[l-(l-x)i/2], (3.35) 
(3.36) 

(3.37) 



Plugging in QCD group factors and doing the integral fl3.35p numerically gives the results 
for a and c in Tables [T] and [III with 



a = 0.55634 Q^/^ 



(3.38) 



The result listed for one-fiavor QED corresponds to setting Ca = = and Cp = tp = 1. 
A similar analysis of the energy stopping distances of fl3.15p gives 



1 



,(<=)«2T V TL 



E 



(3.39) 



where, at leading-log order, 



In 



E 



(3.40) 



as before, with the same value (I3.37P of c. The coefficients a^f \ however, are given by 



Q}/^ {rjiggrnqq^ — mgqmqg) 



'merry TTX/ci 



(3.41) 
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where 



m 



m 



rria 



m 



qq 



(X) 



dx Pq_>g(x) 



Ca + Cax^ + Ca{1-x)^ ^ 
x(l — x) 

(2Cf - Ca) + Cax^ + Ca(1 - xf 



[I _ a;3/2 _ _ ^)3/2] 
1/2- 



x(l — x) 

(2Cf - Ca) + Cax^ + Ca(1 - xf 





a;(l — x) 


Ca + (2Cf 


-Ca)x2 + Ca(1-x)2 




x(l — x) 


Ca + (2Cf 


-Ca)x^ + Ca{1-xY 


x(l — x) 



1/2 



(3.42) 

[_a;3/2 _ (1 _ ^)3/2] ^ (3^43) 



1/2 



1/2 



(3.44) 
. (3.45) 



Results are displayed in Tables HI and HTl 



D. Quark mass thresholds 



In this paper, A^f is the number of effectively massless quark species. Here we comment 
on what counts as a massless quark for the purposes of our analysis. Parametrically, the 
condition to ignore the mass of a high-energy quark in a splitting process q ^ qg or g ^ qq 



IS 



M < ~ (qE) 



1/4 



(3.46) 



So we can use our formulas to describe the quark number stopping time of such a quark, or 
we can compute energy stopping times provided we take N{ in the DGLAP splitting function 
(I3.22P to be the number of quark species that satisfy fl3.46p . All other factors of A^f in this 
paper refer to the number of quark species in the plasma, which should be taken as the 
number satisfying M -C T. The specific results in tables |T] and [TTl are given for the case that 
these two Af values are the same. 



IV. NLLO ANALYSIS OF QUARK NUMBER STOPPING 

Now we extend the analysis of the previous section to next-to-leading logarithmic order 
for the case of the quark number stopping distance. Our goal is to determine the coefficient 
bin (O). 



^ See footnote [TJ For hard bremsstrahlung, the quark mass wiU give a contribution of order JVP / E to 
AE ^ Q\/E. This effect of the mass is neghgible when M <^ Q±. Plugging in the massless result 



Q± (qE)^/* makes this condition ^M^. 
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A. NLLO splitting rates 



The NLLO resuh for sphtting rates was computed in Refs. 11, 13] and corresponds to 
the leading- log result (13.171) . 



dx A7tV2x{1-x)E 
with the formula (I3.18P for /i^^^^ replaced by 



(4.1a) 



2 



4a;(l - x)E 



[-Ca + a + Ce) I {e^'fla-.bc) + " ^ + C,)x^ I 

'e^'fia-^be 



e^'^^a^bc 



X 



+ (a + a-c,)(i 



xY q 



1 — X 



1/2 



where 



(4.1b) 



(4.2) 



e = exp(2-7E + f). 

Numerically, equation (]4.1bl) may be solved self-consistently for fia^bc- Alternatively, one 
may solve it iteratively starting from an initial guess fXa^bc = Q±o ~ {qEY^^ on the right- 
hand side, then using the result as a refined guess, and so forth. This generates an explicit 
expansion of the solution in inverse powers of logarithms. In either case, the final result for 
the rate dV /dx is only correct to NLLO. For the quark-number stopping distance, we only 
care about the q gq splitting rate. 

The explicit expansion of the rate (14. lap in powers of inverse logarithms has the form 



dr 

dx 



I) 



1/2 



A{x) \n^/\nE) 



l\n[\n{KE)] + B{x) 



+ ■ 



(4.3) 



where k is an arbitrary constant with dimensions of inverse energy. B{x) and higher-order 
coefficients in the expansion implicitly depend on k in such a way that the rate is formally 
independent of the choice of k — that is, the dependence on k at any fixed order in the 
expansion is always a higher-order effect. However, in practice k should be chosen to be 
parametrically of order 1/a'^T. (We will make a much more specific choice of k later on.) 

A{x) simply gives the leading-log term, which we will find notationally convenient to 
write as (specializing to g — gq) 



A{x) 



[\ix)] 



1/2 



Anxd 



X] 



with 



X{x) = Qx{l - x) [Ca + (2Cf - Ca)x^ + Ca(1 - xf] 



(4.4) 



(4.5) 



For NLLO, we need better than the leading-log result fl3.25l) for g(A). Complete results 
to leading order in powers of coupling are presented in Refs. jl2]. IT5|.^ For the case A ^ T, 



g(A) ^ a^T^ 



In — 7^ — h A/" In 



TTIt 



(4.6) 



15| also calculated the next-order correction in powers of the coupling which is quite significant 



Ref. 

for any value of as of possible interest to experiment. 
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where 
and 



r = 4exp(l-27E) 



oo ^ 

a+ = V — ln[(A; - 1)!] ~ 0.3860438 



k=l 



k=l 



ln[{k - 1)!] ~ 0.0112168, 



For the case A <^T, the leading-order answer for q is in fact the same as (13.261) 

A2 



g(A) ^ aTm^ In ^ . 



Combining these expressions for q{A) with (14. ip and (14.31) gives 



with 



Bi{x) = iln[2A(x)] -2 



[(2Cf-Ca)x^ lnx + CA(l-x)^ln(l - a:)] 
[Ca + (2Cf-Ca)x2 + Ca(1-x)2] 



B2(X) 



In i. + 4 
In^, 



Q±<^T. 



(4.7) 

(4.8) 
(4.9) 



(4.10) 



(4.11) 

(4.12a) 
(4.12b) 



B. NLLO stopping distance 



Corresponding to the expansion (14. 3 p of the rate, one may now look for a solution to the 
quark number stopping distance equation (13.50 in the form 



^1/2 



1 - 



i ln[ln(/«E)] + 
\n{KE) 



(4.13) 



The internal minus sign in this formula compared to (14.31) arises because the length is related 
to the inverse of the rate. Plugging the ansatz (14.131) and the NLLO rate (14. 3 p into (13. 5p 
determines the same A = aa'^T^^'^ as in the leading-log calculation, with 



a = [ dxA{x) [1 - (1 -x)^/'] . 
Jo 



(4.14) 



It also determines 



1 

B = - dxA{x) {B{x) [l~{l-xf'^] + \{l - xf'Hn{l - x)] 



i ln(a"T/€) + i/? 



(4.15) 
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with 



1 /•! 

f3 = - dxA{x){2[Bi{x) + B2{x)][l-{l-xy/^]+{l-xy^^\n{l-x)}. (4.16) 
a Jo 

As mentioned earlier, the precise choice of k is formally arbitrary, and the effect of 
changing k by a multiplicative factor of order one does not affect the NLLO result except by 
corrections that are yet higher order in inverse logarithms. However, as a practical matter, 
it is useful to have some definite prescription for choosing a sensible value for k. We will 
use the fastest apparent convergence (FAC) prescription, which in this context is to choose 
K such that the NLLO correction in (14.131) vanishes: 

i \n[\n{KE)] - \ ln(a"T/t) + i/? = 0, (4.17) 

and so 

In(KE) = In y ^^^) ■ (4-18) 

This corresponds to fll.3bl) with the identification L = ln{nE) and 

b = e'^. (4.19) 

The combination of equations (I4.12p . (14.161) . and (I4.19P is our final result for the NLLO 
coefficient b in (II. Sp . Evaluating the integrals numerically gives the tabulated results for b 
in Tables [T] and [m 



V. A NUMERICAL TEST OF THE RESULT 

In this paper, we have derived a simple NLLO formula (II. 3p for the quark number stop- 
ping distance of a high-energy particle. This result relies on an expansion in powers of 
inverse logarithms, which parametrically requires ln(ii^/T) 3> 1, and one may wonder just 
how large E has to be before this assumption becomes reasonable. As an alternative, one 
could imagine dropping the assumption that logarithms are large and instead computing 
the stopping distance usin g th e full weak-coupling result for the bremsstrahlung rate, at 
leading order in as 0,1 

llTl . llSl . The disadvantage is that one must then do (somewhat 



complicated) numerics — we know of no simple formula for the stopping distance unless one 
assumes ln{E/T) ^ 1. In this section, we perform such a numerical analysis simply as a 
check that our result (11.30 is indeed correct at sufficiently large E. 

The numerics consist of Monte Carlo evolution of a large sample of quarks with energy 
E, using the full weak-coupling result^ for the bremsstrahlung rate dT/dx to randomly de- 
termine whether each quark loses energy xE in each small time step At. We will focus 
on the case E ^ T and so may (for simplicity) ignore other mechanisms of energy loss, 
such as collisional energy loss. For the same reason, we only consider energy loss and will 



Specifically, we take dTa^bc/dx = {2'r:Y^a^bcl Eva where the 7a^&c are the splitting functions of Refs. 
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18l | and the number Va of spin+color degrees of freedom is 6 for a quark or anti-quark and 16 for a 
gluon. We do not include any final-state Bose enhancement or Fermi blocking factors for the final-state 
particles h and c: In the limit that their energies are large compared to T, no such factors are necessary. 
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ignore processes which can increase the quark energy (such as inverse bremsstrahlung) . 
Also for the sake of simphcity, we will formally focus on the Q±_ -C T case (parametrically 
E ^ T/ctg ln(a~^)), which is the approximation implicitly used in previous numerical cal- 
culations of the bremsstrahlung rate for weak coupling.^ In the numerics, we will call a 
quark "stopped" when its energy drops below T, since our approximations no longer make 
sense then and since a real quark would equilibrate once its energy dropped to be of order 
T. Whether one defines "stopped" as < T or E < T/2 or E < 2T will only affect 
the stopping distance by a relative correction that is parametrically of order ^/tJE (up to 
logarithms) and so is unimportant when E ^ T. 

In each time step, we let each quark bremsstrahlung with probability F At. When 
bremsstrahlung occurs, we randomly choose the momentum fraction of the bremsstrahlung 
gluon according to the probability distribution dT/dx. This is not completely straight- 
forward because the total bremsstrahlung rate F has a small x divergence. Radiation of 
small X gluons is inefficient for energy loss, and so the final result for the stopping distance 
should not be sensitive to this divergence. In our numerics, we simply sidestep such diver- 
gences by restricting the range of x values to ei < a; < 1 — €2, where ei and €2 are small, and 
then we check that our numerical results converge as we make ei and €2 (and At) smaller 
and smaller. 

A comparison of the full numeric result with the NLLO approximation (11. 3p is shown 
in Fig. [2] as a function of E/T for 3-fiavor QCD. The plot shows the size of the relative 
error e = 1 — ^nllo /^numeric that the NLLO approximation makes compared to the numeric 
simulation. The plot extends to absurdly large values of E/T because the purpose of this 
plot is to check our NLLO analysis. We can see that the NLLO result indeed approaches the 
full weak-coupling result at very large E. Moreover, it does not do too badly even at smaller 
values of E/T that are only moderately large: The difference is < 30% for E > 6T. One 
should not take this too seriously, however, since at ~ 6T the various approximations we 
made earlier in this section are in doubt. We do not show any comparison for even smaller 
E/T simply because the NLLO result (11.31) . which implements the FAC prescription, ceases 
to exist: There is no solution to (11.3bp for L when E/T < a'^e/b, which in the present case 
is E/T < 5.74. At exactly E/T = a^e/b, the logarithm L becomes 1, and so the idea of a 
large-logarithm approximation has certainly broken down then. 

As a further check of our NLLO calculation, we plot eln^(i?/T) versus y = 1/ \n(E/T) 
in Fig. [3l If the NNLO contribution was incorrect, the result would diverge a.s y 0. 

Finally, just to give some idea of the distribution of stopping lengths in the high-energy 
limit, we plot the distribution from our highest-energy numerical simulation in Fig. HI For 
some quantities of interest in the subject of energy loss (such as the medium effect on energy 
loss from passage through a thin medium), there are slowly falling tails to distributions which 
cause average values of those quantities to be unrepresentative of typical events. That is not 
the case here for the stopping length: the average coincides with the location of the peak of 
the distribution. 



^ Specifically, numerical calculations such as Refs. 0, [l^l have used the formula A{q±) = 'nv^Tlq\[q\ + 
m^) in the notation of Ref. [tsj. ( A{q \ ) corresponds to {27: / gYdTc\l<Pqi_ in this paper's notation and 
T C{q±) in the notation of Ref. [l9|.) This formula is correct only when <C T. As noted in Ref. 11|, 
for 3-flavor QCD, use of this formula in the Q± ^ T case should produce an error of at most 15%. 
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FIG. 2: The relative difference e = 1— tNLLo/^numeric between (i) tlie NLLO quark number stopping 
distance (|1.3p and (ii) tlie quark number stopping distance tnumeric computed tlirougli simulation 
using the full weak-coupling result for the bremsstrahlung rate. The slight scatter of points is due 
to the statistics of simulating only a finite number of particles (in this case, 10^ particles). 
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FIG. 3: The results of Fig. [2] plotted with axes convenient for checking the correctness of the 
NLLO result at very large E. 



VI. RUNNING COUPLING 



Our result (11.31) formally assumes that a\n{E/T) <^ 1 and so ignores, for instance, 
running of the coupling a. One can plausibly accommodate running in the bremsstrahlung 



rate, and so the stopping distance, by the prescriptions reviewed in Ref. 12|.^ The in 



the denominator of (11.31) comes from (i) the overall factor of a in (I3.17P associated with 
the emission vertex for the bremsstrahlung gluon, and (ii) the factor /i^ in (13.171) . which 



is proportional to ^/q oc yc?. As reviewed in Ref. [12], the emission vertex factor should 



See also the earlier discussion in Refs. 
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FIG. 4: Simulation results for the probability distribution for quarks to lose their energy through 
bremsstrahlung in a distance i, as function of i/ia.vg, where £avg is the average. The particular 
results shown are for E/T = 10^, for which £avg agrees with the NLLO stopping length (|1.3p to 
within 4%. 



plausibly be evaluated at the scale Q± ~ /i ~ and the running of a in g replaces 

g(A) oc — > g(A) oc a(A) a(mD). (6.1) 
Since we needed g at A ~ Q^, the effect of both these replacements is to replace (II. 3p by 



1 / E 



''Stop,q 



aa3/2(Q^)ai/2(mo)T V TL 



(6.2) 



with ~ {qEY/\ The origin of the in (]1.3bl) is the difference between the scales 
mD and T in the arguments of the logarithms of (14.61) . and so that a'^ should simply be 
[a(mD)]'^. Throughout this discussion, we do not distinguish between a(T) and a{mD) since 
their difference is small when a is small. 

The result fl6.ip only accounts for aln{E/T) corrections associated with running of the 
coupling. We are not sure whether there might be other corrections of comparable size, 
which is why our introduction is only so bold as to claim a definite result in the formal limit 
aln{E/T) < 1. 
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This is a slight simplification. This nice compact formulation of the prescription only works if there are 
no quarks masses between mu and A. 
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